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Abstract 

The research airdata system of an instrumented F-104 
aircraft has been calibrated to measure winds aloft in sup- 
port of the space shuttle wind measurement investigation at 
the National Aeronautics and Space Administration Ames 
Research Center’s Dryden Flight Research Facility. For 
this investigation, wind measurement accuracies compara- 
ble to those obtained from Jimsphere balloons were de- 
sired. This required an airdata calibration more accurate 
than needed for most aircraft research programs. The F-104 
aircraft was equipped with a research pitot-static noseboom 
with integral angle-of-attack and flank angle-of-attack vanes 
and a ring-laser-gyro inertial reference unit Tower fly- 
bys and radar acceleration-decelerations were used to cal- 
ibrate Mach number and total temperature. Angle of at- 
tack and angle of sideslip were calibrated with a trajec- 
tory reconstruction technique using a multiple-state lin- 
ear Kalman filter. The F-104 aircraft and instrumentation 
configuration, flight test maneuvers, data corrections, cal- 
ibration techniques, and resulting calibrations and data re- 
peatability are presented. Recommendations for future air- 
data systems on aircraft used to measure winds aloft are 
also given. 


Nomenclature 

Values are given in the English system of units. 
Acronyms 

AICS airborne instrumentation compu- 

ter system 

c.g. center of gravity 

IRU inertial reference unit 


•Aerospace Engineer. 
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LKF 

linear Kalman filter 

PCM 

pulse code modulation 

Symbols 

A 

acceleration, g 

At 

normal acceleration at noseboom flow 
vanes, ft/sec 2 

9 o 

acceleration due to gravity, 32.1740 ft/ 

hp 

pressure altitude above mean sea level, 

k 

temperature recovery factor 

M 

Mach number 

P. 

static pressure, lb/ft 2 

Ptot 

total pressure, lb/ft 2 


P 

V 

Q 

Q 

9 

R 

r 

r 

T 

u 

V 
v 

W 

w 

X 


roll rate, rad/sec 
roll acceleration, rad/sec 2 
calculation variable 
pitch rate, rad/sec 
pitch acceleration, rad/sec 2 
ideal gas constant, 53.35 
yaw rate, rad/sec 
yaw acceleration, rad/sec 2 
temperature, °R 
airspeed along X axis, ft/sec 
aircraft velocity, fi/sec 
airspeed along Y axis, ft/sec 
wind velocity, ft/sec 
airspeed along Z axis, ft/sec 

distance from aircraft c.g. along aircraft 
X axis, ft 

distance from flow vane along aircraft 
X axis, ft 



V 

m 

v 

z 

z 


z 


distance from aircraft c.g. along aircraft 

Y axis, ft 

distance from flow vane along aircraft 

Y axis, ft 

geometric altitude above mean sea level, ft 

distance from aircraft c.g. along aircraft 
Z axis, ft 

distance from flow vane along aircraft 
Z axis, ft 


Greek 

a 

P 

r 

7 

6 

9' 

O'g 

4> 

4>' 


angle of attack, deg 
angle of sideslip, deg 
Euler rotation transformation matrix 
ratio of specific heats for air, 1.4 
pitch angle, deg 

noseboom misalignment in pitch, deg 

noseboom deflection in pitch due to 
elevated g, deg /g 

roll angle, deg 

noseboom misalignment in roll, deg 

heading angle, deg 

noseboom misalignment in yaw, deg 


Superscripts 


T transpose 

* reference value calculated by LKF 


Subscripts 

1 

2 

a 

bias 

c 

D 

E 

F 

fb 

I 

i 

m 

N 

n 

t 

tot 


corrected for misalignment 

corrected for misalignment and 
angular rates 

in aircraft axes 

Y -intercept of angle error as a function 
of corrected angle 

all analytic corrections made 

down 

east 

flank 

at fly-by tower reference altitude 

inertial 

indicated 

derived by meteorlogical analysis 
north 

in noseboom axes 

true 

total 


z Z axis accelerometer 

a angle-of-attack vane 

p flank angle-of-attack vane 

oo free stream 

Introduction 

The most common method of obtaining atmospheric wind 
profiles involves tracking a rising balloon (Jimsphere or 
Rawinsonde) which moves with the winds. Balloon meth- 
ods are adequate for many applications but have some limi- 
tations. For example, it is impossible to control a balloon’s 
flight path once launched. In addition, the typical rise rate is 
approximately 15 ft/sec (5 m/sec), so approximately an hour 
is required to obtain an altitude of 60,000 ft (18,000 m). For 
certain applications, such as the space shuttle program, it is 
desirable to obtain wind profiles quickly, in approximately 
10 to 15 min. 

The preprogrammed launch trajectory for the space shut- 
tle is based partially on an expected wind profile for the time 
of year. If the winds on the day of launch are significantly 
different from the expected winds, certain shuttle structural 
load limits may be exceeded. Currently, winds are mea- 
sured on the day of launch by a series of Jimsphere bal- 
loons. The last balloon for loads assessment is launched 2 
hr before shuttle launch. After the balloon reaches 60,000 ft 
(18,000 m), the wind data are used in trajectory simulation 
and loads prediction programs. A wind persistence factor is 
added to the calculations to account for wind changes over 
time and the possibility that the balloon has blown away 
from the launch path. Studies have shown that these changes 
increase significantly when delays are longer than 2 hr or 
spatial separations are greater than 12 mi (20 km). 1-3 De- 
pending on the results of the loads predictions, a go-no-go 
recommendation for launch is made. 

A reduction in the uncertainties in the prelaunch wind 
load assessment caused by temporal and spatial variabil- 
ities was desired. The National Aeronautics and Space 
Administration (NASA) Johnson Space Center requested 
NASA Ames Research Center’s Dryden Flight Research Fa- 
cility to perform a flight experiment to determine the fea- 
sibility of using an instrumented high-performance aircraft 
to measure wind profiles. For this technique to be appli- 
cable to the space shuttle program, the following guide- 
lines were suggested: (1) obtaining a profile to 60,000 ft 
(18,000 m) in 10 min, (2) an aircraft profile area within a 
10 mi (16 km) radius circle, and (3) measurement accu- 
racy comparable to the Jimsphere system. Winds are cal- 
culated as the vector difference between aircraft inertial ve- 
locity, measured by an onboard inertial reference unit (IRU) 
or ground based radar, and aircraft airspeed, measured by a 
calibrated airdata system. 

The quality of the wind profiles measured by Jimsphere 
balloons can be shown by the tracking precision for an in- 
dividual balloon and the repeatability of winds from sepa- 
rate balloons. Two separate radar systems’ trackings show 
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a root-mean-square (rms) error for the Jimsphere balloon 
of 1.6 ft/sec (0.5 m/sec) for wind velocities averaged over 
50 m intervals. 1 To match this precision, an aircraft traveling 
at Mach 2 would need to measure its Mach number within 
0.0016 and its flow angles to within 0.05°. Studies show that 
the rms repeatability of wind profiles from two Jimsphere 
balloons is approximately 7 to 10 ft/sec (2 to 3 m/sec) if the 
time separation is less than 1 hr or the distance separation 
is less than 12 mi (20 km). 1-3 This is comparable to 0.007 
to 0.010 in Mach number and 0.2° to 0.3° in flow angles 
at Mach 2. (Actual aircraft wind measurement repeatability 
depends on the accuracy of several sensors as well as on the 
atmospheric variability between wind measurements.) Be- 
cause of these tight tolerances on Mach number and flow 
angles, an airdata calibration more accurate than required 
for most aircraft research programs is necessary. 

This paper covers the methods and results of the airdata 
calibration of an F- 104 aircraft used to measure winds aloft 
in support of the F- 104 shuttle wind measurement investiga- 
tion. Wind profile data were gathered at altitudes from 3,000 
to 67,000 ft (900 to 20,400 m), Mach numbers from 0.78 to 
2, angles of attack from 0° to 12°, and angles of sideslip from 
-3° to 3°. The F-104 aircraft and instrumentation config- 
uration, flight test maneuvers, data corrections, calibration 
techniques, and resulting calibrations and data repeatabil- 
ities for Mach number, total temperature, and flow angles 
are presented in this report. Recommendations for future 
airdata systems on aircraft used to measure winds aloft are 
also given. 

Background 


ment collar. Indicated static pressure was measured by the 
forward set of static pressure orifices; the pilot’s instrumen- 
tation and the aircraft airdata computer used the aft set (Fig. 
2(b)). Total pressure and static pressure were each measured 
with a separate high-resolution, high-accuracy absolute- 
pressure transducer. Angle-of-attack and flank angle-of- 
attack values from the vanes were measured using poten- 
tiometers. Total temperature was measured by a nonde- 
iced, open element sensor (Fig. 3). Research parameters 
were digitally encoded using pulse code modulation (PCM) 
and were both recorded onboard and telemetered to ground 
based recorders. 

A commercial airline aircraft type ring-laser-gyro ERU 
was used to give inertial velocity, acceleration, attitude, and 
angular velocity. 11 An airborne instmmentatfflfi computer 
system (AICS) interfaced the IRU with the PCM system 
and supplied the IRU with inputs from the airdata system. 12 
These inputs are needed to stabilize the IRU integration 
loop. Figure 4 shows the IRU, AICS, PCM system, and sig- 
nal conditioning hardware installed in the aircraft. 

Since the IRU was designed for commercial airline air- 
craft, its parameters, especially normal acceleration, were 
heavily low pass filtered. Because high rate maneuvers were 
used for the flow angle calibrations, it was decided to use a 
separate accelerometer set located behind the cockpit The 
aircraft also had equipment which allowed trajectory guid- 
ance information to be telemetered from the ground and 
displayed to the pilot. 13 This guidance method is known 
as uplink. 

Procedure 


Aircraft have been used to obtain wind measurements, 
but most results have been obtained during steady-level 
flight. 4-6 In steady-level flight, many simplifying assump- 
tions are usually made. These assumptions are no longer 
valid in high-speed, climbing, cm- descending flight, as used 
to minimize time and spatial separation for a shuttle appli- 
cation. The wind measurement methods for an aircraft in 
a nonsteady flight environment have been documented in 
Ref. 7. 

Traditional airdata calibration methods are well 
known, 8 - 9 but a higher accuracy than is typically achieved 
with these methods was desired. The methods used to 
achieve these accuracies are described in this report. 


Mach Number Calibration 

Mach number was calibrated using tower fly-by and radar 
acceleration-deceleration methods. For the tower fly-by 
method, shown in Fig. 5, the aircraft is flown at a steady 
airspeed and altitude near the fly-by tower. 9 The aircraft is 
sighted from the tower through an eyepiece and grid and the 
aircraft’s true geometric altitude ( Z ) is determined by geom- 
etry. The ffee-stream static pressure at the aircraft’s altitude, 
( P'J is calculated using the formula 


P.„ = 


e 




(i) 


Description of Aircraft and Instrumentation 

The aircraft used for the wind-measurement flight tests 
was an instrumented F-104 (Fig. 1). A standard National 
Advisory Committee for Aeronautics (NACA) research 
noseboom (Fig. 2(a)) was used to measure static and to- 
tal pressures and angles of attack and sideslip. 10 Figure 2(a) 
also shows the size and location of the noseboom attach- 


where P, /t , T/ h , and Z/t, are the static pressure, am- 
bient air temperature, and geometric altitude at the fly- 
by tower reference altitude, respectively. This for- 
mula is the integral form of the hydrostatic equation. 14 
The total pressure measured by the noseboom (P totj ) 
is assumed to be correct, 8 and the free-stream Mach 
number (Moo) for subsonic flight is calculated using 
the expression 
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where 


Mac 



( 2 ) 


Q = 1.839371 



Indicated Mach number (Mi) is calculated using equa- 
tion (2) with indicated static pressure ( P *) substituted for 
free-stream static pressure. Passes by the tower are flown at 
various subsonic Mach numbers, and the Mach correction 

AM = Moo — Mi (3) 

is plotted against indicated Mach number to give the Mach 
calibration. The calibration was assumed to be a function 
of indicated Mach number only, since high-lift-coefficient 
flight was not part of the investigation. The tower fly-by 
technique is the most accurate method, but only subsonic 
data can be taken, and only a small number of calibration 
points can be flown on one flight. 

In the radar acceleration-deceleration technique (Fig. 6), 
as in the tower fly-by technique, free-stream static pressure 
is calculated for the aircraft. For a calibration run, the air- 
craft flies with wings level on a constant heading at a con- 
stant geometric altitude. The aircraft’s radar altitude is up- 
linked to the aircraft, and the difference between the desired 
and actual geometric altitude is displayed in the cockpit to 
aid the pilot in flying a constant geometric altitude. The 
aircraft begins the run at a low airspeed, and accelerates at 
approximately 3 kn/sec to a peak Mach number. The pilot 
then begins to decelerate at approximately 3 kn/sec back to 
the original low airspeed. The entire maneuver is executed 
at radar elevation angles above 10° to minimize radar re- 
fraction errors and below 80° to avoid high radar antennae 
slew rates. 15 Time coded radar data is processed to give geo- 
metric altitude. Weather balloon and atmospheric chart data 
are analyzed to obtain a table of geometric altitude minus 
pressure altitude ( Z — hp) m as a function of geometric alti- 
tude (Z). The pressure altitude error caused by lateral atmo- 
spheric pressure gradient (A hp m ) with lateral distance and 
direction from the radar site is also tabulated as a function of 
geometric altitude. Applying these tables to the time history 
of the aircraft’s radar altitude yields the true pressure alti- 
tude of the aircraft, which is then converted into free-stream 
static pressure. 16 As with the tower fly-by method, indi- 
cated total pressure is assumed to be true in subsonic flight, 
and only normal shock losses are assumed in supersonic 
flight 8 Equation (2) is then used to calculate free-stream 
Mach number subsonically, (when < 1.89293, 

the value at M «, = 1.0). 


is used supersonically, (when PtotJP» x > 1.89293). This 
equation is a Thylor series expansion of the Raleigh pitot- 
static equation. 16 Indicated Mach number is calculated by 
using equations (2) and (4) with indicated static pressure re- 
placing free-stream static pressure. As with the tower fly-by 
data, equation (3) is used to obtain the Mach correction. The 
result is plotted against indicated Mach number. The radar 
acceleration-deceleration method is useful because: (1) it 
allows supersonic Mach calibration, (2) it is independent of 
the tower fly-by method, (3) the entire Mach calibration can 
be obtained in a few minutes of flight time, and (4) the ma- 
neuvers can be conducted at several altitudes to check for 
altitude effects on calibration. 

The radar acceleration-deceleration Mach corrections are 
inconsistent with those derived by the tower fly-by method. 
This difference in Mach correction is a pure bias in altitude. 
Compared to the tower fly-by method, the pressure altitude 
derived from radar and weather analysis is too low. The lo- 
cus of points, obtained on days when the atmosphere was 
relatively calm, has a generally increasing magnitude of al- 
titude bias with increasing altitude. The bias ranges from 0 
error near the ground to as much as —264 ft at a pressure 
altitude of 34,000 ft. The cause of this altitude bias is being 
investigated, but is thought to be a systematic error in the 
measurement or analysis of the weather or radar data. 

The altitude bias error in the radar acceleration- 
deceleration data was eliminated by “bootstrapping.” For 
this technique, the indicated Mach number for a single, sub- 
sonic point at the beginning of the acceleration is calcu- 
lated. The pressure altitude is calculated using the Mach 
calibration derived from a tower fly-by. The difference be- 
tween this pressure altitude and that calculated from radar 
and the (Z - hp) m and A hp m tables from the weather 
analysis is the altitude bias. The altitude bias correction 
is then applied to the entire maneuver. The correction is 
valid as long as the maneuver is flown at a nearly con- 
stant altitude. The resulting calibration agrees well with 
the tower fly-by calibration over the entire subsonic Mach 
range, even though only one point of the radar acceleration- 
deceleration calibration was forced equal to the tower fly- 
by calibration. At a given Mach number the aircraft al- 
titude is correlated to the aircraft angle of attack. Since 
the static pressure measured by a standard NACA airdata 
probe is insensitive to angles of attack up to 20°, 10 the 
aircraft’s Mach number calibration should not depend on 
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angle of attack, or thereby altitude. Therefore, the Mach 
number calibration should not be a function of the altitude 
bias correction. 

Total Temperature Calibration 

The ambient air temperature at the aircraft altitude (T^) 
may be calculated from the measured total temperature 
CTtot), and free-stream Mach number by 



where k is the recovery factor of the total temperature sen- 
sor. The value of k is usually a constant parameter and 
is dependent on the sensor geometry. The recovery factor 
can be calculated using tower fly-by and radar acceleration- 
deceleration data. If total temperature is plotted as a func- 
tion of free-stream Mach number squared, the slope is 
Too k, and the y-intercept is the free-stream temper- 
ature. A least-squares curve fit was used to determine 
these quantities (Fig. 7). The value of the recovery fac- 
tor for the aircraft is taken to be the average from all the 
maneuvers analyzed. 

Flow Angle Corrections 

The vanes on the noseboom measure the angles between 
the local velocity vector and the noseboom axes. Figure 8 
shows the steps taken to correct these angles to true angles of 
attack and sideslip: (1) noseboom misalignment correction, 
(2) aircraft angular-rate correction, (3) noseboom bending 
correction, (4) aerodynamic flow angle calibration correc- 
tion, and (5) transformation of flank angle of attack to angle 
of sideslip. The first, second, third, and fifth steps are analyt- 
ical corrections. The fourth is determined through trajectory 
reconstruction techniques. Figure 9 shows the difference be- 
tween flank angle of attack and angle of sideslip. 

Noseboom Misalignment Corrections. The first cor- 
rection to the flow angles (Fig. 8) is for noseboom mis- 
alignment. The noseboom on an aircraft may be imperfectly 
aligned with the aircraft axes. To transform the angles from 
the noseboom axes system (X n , Y n , and Z n ) to the aircraft 
axes system (X a , Y a , and Z a ) (Fig. 10), the three compo- 
nents of the true airspeed in the noseboom axes system (u n , 
v n , and w n ) must be calculated. The free-stream velocity (in 
ft/sec) is calculated by 


V«, = M 00 \/'yRgoT 00 (6) 


The three components of airspeed are then evaluated 


Voo 

\/l + tan 2 (a,) + tan 2 (a/.;) 


( 7 ) 


v„ = u n tan(a Fj ) 
w n = tin tan (a,) 


where a, is the indicated angle of attack and a Pi is the indi- 
cated flank angle of attack. The angular offset of the nose- 
boom in roll (<p’), pitch (O'), and yaw (ip') is measured from 
the aircraft axes to the noseboom axes (Fig. 10). An Euler 
rotation through these angles transforms the noseboom axes 
(u„, v„, and id,) components into the aircraft axes (u„, v„, 
and w a ) 


where 




x 


cos(V>') — sin(V>') 0 

sin (VO co s(rp') 0 
0 0 1 

cos(0') 0 sin(0') 

0 1 0 
-sin(<?') 0 cos (O') 


x 


1 0 0 

0 cos(<p') -sin (<p') 

0 sin(<£') cos (<p') 


( 8 ) 


The indicated angle of attack and flank angle of attack can 
be expressed in the aircraft axes system by 


ai=tan -1 (*) 

, ) \ < 9 > 

a*=tan-‘ (*) 

The axes system transformation was necessary because of 
the high accuracy of the angle-of-attack and -sideslip mea- 
surements needed to accurately measure winds. At an an- 
gle of attack of 20° and an angle of sideslip of 5°, on this 
aircraft, neglecting the noseboom misalignment correction 
would result in a 0.5°-enor in flow angle. 

To measure the boom offset angles for roll, pitch, and yaw 
a clinometer and transit are used. The roll and pitch attitude 
of the aircraft are measured with the clinometer at reference 
points on the aircraft in the hangar. The clinometer is then 
used to measure the pitch and roll angles of the noseboom, 
taking into account the 0.33°-taper in the support shaft of 
the angle-of-attack vane. The differences are the angular 
offset for roll and pitch respectively. The yaw offset angle 
is determined by sighting two hard points on the aircraft and 
two points on the noseboom using the transit 

The noseboom is fastened to the aircraft radome, which is 
removed periodically for maintenance. After removing and 
reinstalling the radome, the boom misalignment measure- 
ment is repeated. The values for the offset angles for two 
flight periods are given in Table 1. These quantities vary by 
a significant amount 

Aircraft Angular-Rate Corrections. The second cor- 
rection to the flow angles (Fig. 8) is for angular rates of the 
aircraft. As the aircraft rolls, pitches, and yaws, the motion 
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of the noseboom about the center of gravity (c.g.) induces 
additional airspeed components, deflecting the flow vanes. 
For angle of attack, the pitching and rolling effects of the 
aircraft were accounted for by 

. _ : _-l ((qXa-PV a )COS(cn)\ /im 

a 2 = ai + sin I — I ( 10) 

where p, q, and r are the aircraft roll, pitch, and yaw rates, 
respectively, and x a and y a are the longitudinal and lateral 
distances from the aircraft c.g. to the angle-of-atlack vane. 
These distances are tabulated in Thble 2. For flank angle of 
attack, rolling and yawing effects are accounted for by 


cap e — oiF x + sin 


_1 { (pzp - rx p) COS ( Off, ) \ 

V ) 


(11) 


where zp and xp are the normal and longitudinal distances 
from the aircraft c.g. to the flank angle-of-attack vane. 
These are also given in Table 2. Figure 11 shows the 
angular-rate corrections. Equations (10) and (11) are not 
limited by small angle approximations and are applicable to 
high angle-of-attack and high angle-of-sideslip calculations. 

These angular rates, like other parameters measured by 
the IRU, have time delays and need to be skewed the appro- 
priate amount in time. The time delays for this configuration 
are shown in Table 3. 


Noseboom Bending Correction. The third flow an- 
gle correction (Fig. 8) is for noseboom bending. During 
elevated-g maneuvering, the noseboom deflects under the 
increased load. The noseboom deflection in pitch caused by 
elevated g (ff g ) was measured experimentally for this aircraft 
as -0.064%. No correction is made for noseboom bending 
in the yaw plane, since the aircraft is not significantly loaded 
in that plane. 


The aircraft’s accelerometer package is located behind the 
cockpit in the electrical bay, so acceleration data must be 
translated to the noseboom before the noseboom deflection 
can be calculated. Thble 2 presents the distances x„y„ and 
z„ from the noseboom to the vertical accelerometer (/I*). 
The normal acceleration at the noseboom (At) is 


At = —A, + 


(pr-g)x, + (gr+p)p, -(p 2 + q 2 )z M 


90 


( 12 ) 


The roll, pitch, and yaw accelerations (p,q, and f) were nu- 
merically differentiated from p, q, and r using a five point 
weighted least-squares sliding window. 17 

Angle of attack corrected for noseboom bending is 

a c = ctj — O'giAt — cos(O) cos(<f>)) (13) 

where the term involving the cosines subtracts the normal 
acceleration caused by gravity. 


Aerodynamic Flow Angle Calibrations Corrections. 
The fourth correction to the flow angles (Fig. 8) is for aero- 
dynamic effects that the aircraft and the noseboom induce on 


the local velocity vector. These effects have been calibrated 
through flight test data using a trajectory reconstruction al- 
gorithm. The corrections for angle of attack will be called 
upwash and for flank angle of attack the corrections will be 
called sidewash. 

For subsonic flight, upwash was identified using the fol- 
lowing flight test technique. Aided by the uplink, the pi- 
lot held the aircraft at a constant altitude, Mach number, 
and power setting for several seconds, then swept through a 
range of angle of attack, followed by several seconds of sta- 
bilized data (Fig. 12(a)). Sidewash was calibrated in a simi- 
lar fashion with sweeps in flank angle of attack (Fig. 12(b)). 
These maneuvers were conducted at true Mach numbers 
ranging from 0.78 to 0.92. 

Since no flow disturbances from the aircraft can propa- 
gate forward to the flow vanes in supersonic flow, upwash 
and sidewash will not be a function of angle of attack or 
flank angle of attack. Shock waves on the noseboom af- 
fect the indicated flow angles however, so upwash and side- 
wash are identified as a function of indicated Mach num- 
ber, using the radar acceleration-deceleration maneuver de- 
scribed earlier. 

The true angle of attack and flank angle of attack are 
calculated by a trajectory reconstruction algorithm with a 
multiple-state linear Kalman filter (LKF). 18 The LKF blends 
data from the aircraft’s accelerometers, IRU, airdata system, 
radar tracking, and weather analysis to give the minimum 
variance estimate of the aircraft’s trajectory. The observa- 
tions and the dynamics equations are selectively weighted 
using a matrix determined by physical intuition about the 
system. Tables 4 and 5 give the observation weighting ma- 
trixes used for the upwash and sidewash calibrations. Ta- 
bles 6 and 7 show the weightings used for the dynamics 
equations. The LKF algorithm consists of a prediction and 
a correction step. The prediction step extrapolates the mea- 
sured data to the next time point using the dynamics equa- 
tions. The correction step adjusts the extrapolated state us- 
ing measured data at that next time point to give the mini- 
mum variance estimate. 18 

The components of the inertial velocity and winds aloft, 
calculated by the LKF, were added in vector form to give 
estimates of true angle of attack and flank angle of attack 
for each time point in the calibration maneuver, that is at 
each data frame. True angle of attack and true flank angle 
of attack are determined from the LKF reference states as 


a t = tan -1 



(14) 


_i v 

atp, = tan r 
u 


( 15 ) 


where 


u 


f 

' v N ' 


' W N ■ 

> 

V 

= r T . 


Ve 

— 

W E 


w 



V d 


W D 
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r T = 


1 o 

0 cos(<f>) 
0 -sin(<f>) 


x 


cos(0) 

0 

sin(0) 


0 

1 

0 


0 

sin(<£) 

cos(<f>) 

- sin(0) 
0 

cos (0) 


x 


cos (VO sin (VO 0 
-sin(^) cos(V0 0 
0 0 1 


(16) 


There are several things to note in these equations. First, 
true flank angle of attack is used since it is independent of 
angle of attack and more directly comparable to corrected 
flank angle of attack than true angle of sideslip is. Second, 
the vertical component of the wind measurement is set to 
zero since it is difficult to measure using weather balloons. 
Lastly, r T is used to transform the velocity components 
from earth axes to aircraft body axes. 

The errors in angle of attack and flank angle of attack are 


A a = a c — at 

(17) 

Aa? = ap c — ap t 

(18) 


These errors are plotted as functions of corrected angle of 
attack, corrected flank angle of attack, and indicated Mach 
number in order to identify systematic trends that can be 
used for a calibration. Figure 13 shows typical subsonic cal- 
ibration maneuvers at a free-stream Mach number of 0.92. 
Once identified, a line is fit to the data using least-squares 
regression. The slope of the error in angle of attack as a 
function of corrected angle of attack curve ) is defined 

as the upwash factor, and the slope of the error in flank an- 
gle of attack as a function of corrected flank angle of attack 

curve is the sidewash factor. The y-intercepts of 

the curves are the angle-of-attack and flank angle-of-attack 
biases (Qbm and a^). The upwash and sidewash factors 
are a function of indicated Mach number as is the angle-of- 
attack bias. 


No trend of flank angle-of-attack bias with indicated 
Mach number was found in this calibration. The high sen- 
sitivity of flank angle of attack to small changes in the wind 
estimate makes this task impossible without many more data 
flights to get a statistically significant set of data points. This 
is a limitation of the method since it relies on an accurate 
measurement of the winds by weather balloons. 

Subsonically, the true angle of attack and true flank angle 
of attack are calculated by 


a t 


a e — 


6 A a 

~z a c + Qffeiss 

oa c 


&F, = a Fc - 


6 A ap 

— ap + a Fhim 

Sa Fc 


(19) 

( 20 ) 


for A/,' < 1.0. 


As noted previously, in supersonic flight flow angle er- 
rors are a function of Mach number only. The trend with 
Mach number is continuous and repeatable outside of the 
transonic region, where local shock waves interfere with the 
flow vanes. 19 True angle of attack and true flank angle of at- 
tack for supersonic flight are 


a ( = a c -Aa 

(21) 

a p t = a p c - A ap 

(22) 


for M,> 1.0. 

The calibration coefficients can be used as long as the ex- 
ternal configuration of the aircraft is not significantly altered 
and the sensors have not degraded. It is important to have 
preflight inspections of the noseboom, total temperature sen- 
sor, IRU, and supporting electronics to assure quality data 
will be gathered. 

Transformation of Flank Angle of Attack to Angle of 
Sideslip. The fifth and last correction to the flow angles 
(Fig. 8) is the transformation of flank angle of attack to angle 
of sidelip, given by 

0t = tan -1 [tan(af,) cos(ar t )] (23) 

The distinction between these two angles is shown in 
Fig. 9. 


Results and Discussions 

Mach Number Calibration 

To calibrate Mach number for this aircraft, 22 tower fly- 
by points and 11 radar acceleration-deceleration maneuvers 
were obtained and analyzed. Figure 14 shows the Mach 
number calibration correction as a function of indicated 
Mach number for all the calibration maneuvers. All of the 
radar acceleration-deceleration data were corrected using 
the bootstrap technique. As Fig. 14 shows, the scatter is ap- 
proximately ±0.003 in the subsonic data, and approximately 
±0.005 in the supersonic data. A curve faired through the 
data was used as the Mach number calibration, shown in 
Fig. 14 as a solid line. 

Total Temperature Calibration 

Total temperature was calibrated by extracting the re- 
covery factor from tower fly-by and radar acceleration- 
deceleration data. Table 8 shows the individual recovery 
factors with magnitudes varying from 0.977 to 0.995, with 
an average of 0.986. Some of the ±0.009 scatter is caused 
by variations in ambient air temperature during the calibra- 
tion maneuvers, especially for the tower fly-by maneuvers 
which were conducted over more than an hour of flight time. 
Wind-tunnel data for this type of temperature sensor yields 
recovery factors between 0.992 and 0.999. 20 

Flow Angle Calibrations 

The flow angles for this aircraft were corrected using the 
five steps shown in Fig. 8. One of these steps was correcting 
the flow angles for aerodynamic upwash and sidewash. 
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Subsonic Upwash Calibration. Figure 15 shows the 
subsonic upwash factor as a function of indicated Mach 
number. TTie increase of the upwash factor with increas- 
ing Mach number is thought to be caused by a compression 
effect on the angle-of-attack vane, caused by the close prox- 
imity to the flank angle-of-attack vane support shaft Dual 
flank angle-of-attack vanes, one opposite the original vane, 
would make the compression field symmetrical about the 
X n - Y n plane, and negate the change in upwash factor with 
indicated Mach number. The dependance of the angle-of- 
attack bias on indicated Mach number is shown in Fig. 16, 
along with its calibration fairing. Comparing the fully cali- 
brated subsonic angle of attack with the LKF estimate yields 
a residual bound of approximately ±0.2°. 

Subsonic Sidewash Calibration. The subsonic side- 
wash factor plotted as a function of indicated Mach num- 
ber is shown in Fig. 17. The increase in sidewash factor 
with increasing Mach number above = 0.84 is thought 
to be caused by a compression effect on the flank angle-of- 
attack vane caused by the noseboom attachment collar lo- 
cated immediately behind the vane (Fig. 2(a)). The flank 
angle-of-attack bias had too much scatter to calibrate be- 
cause of its high sensitivity to errors in the meteorological 
winds estimate. For these calibrations, the flank angle-of- 
attack bias was assumed to be zero. Comparing the fully 
calibrated subsonic flank angle of attack with its LKF ref- 
erence estimate yields a residual whose variations about the 
mean are bounded by ±0.15°, meaning sidewash factor has 
been calibrated to within ±0. 15°. 

Supersonic Upwash Calibration. Figure 1 8 shows the 
supersonic angle-of-attack error as a function of indicated 
Mach number for several flights with its faired calibration 
curve. The calibration curve shows a sharp drop in error at 
Mi = 1.20, which is thought to be caused by a local shock 
wave crossing the angle-of-attack vane. At higher indicated 
Mach numbers, angle-of-attack error is constant As Fig. 18 
shows, the flight-to-flight scatter in the supersonic angle-of- 
attack calibration is approximately ±0.3°. 

Supersonic Sidewash Calibration. Figure 19 shows 
the supersonic flank angle-of-attack error as a function of in- 
dicated Mach number for several flights along with its faired 
calibration curve. The decrease in the change in flank an- 
gle of attack above A/, = 1.250 is thought to be caused by 
the single angle-of-attack vane. The oblique shock wave off 
of the angle-of-attack vane and support shaft increases the 
static pressure on the port side of the noseboom, making the 
flank angle-of-attack vane read too low. This effect becomes 
more pronounced as indicated Mach number increases be- 
cause of the increasing shock strength, as the curve’s neg- 
ative slope shows. Dual angle-of-attack vanes would make 
the pressure distribution around the noseboom symmetric 
about the X n - Z n plane, negating the Mach effect on flank 
angle of attack. As Fig. 19 shows, the flight-to-flight scatter 
in the supersonic flank angle-of-attack calibration is approx- 
imately ±0.25°. 


Effect of Airdata Calibration on Wind 
Measurement Quality 

Assuming negligible errors in the inertial velocities cal- 
culated by the IRU or ground based radar, the calibration re- 
peatabilities attained should result in measured wind speed 
repeatabilities of 3 to 10 ft/sec (1 to 3 m/sec). The airdata 
repeatabilities, in terms of the calibration error bounds, are 
similar to the rms repeatability as the radar tracked Jim- 
sphere balloon. 1 ~ 3 The wind speed repeatability of 3 ft/sec 
is calculated assuming a Mach number of 0.9, and the 
10 ft/sec repeatability assumes a Mach number of 2. Since 
the flow angles are used to resolve the airspeed into wind 
components, a given error in a flow angle will cause an in- 
creasing wind error as aircraft speed increases. For this air- 
craft high speeds, however, are necessary to obtain high al- 
titudes and to rapidly survey the altitude desired. 

Recommendations 

Several recommendations from this experiment may be 
helpful in the design and operation of future airdata sys- 
tems for aircraft used to measure winds aloft To minimize 
the effects of Mach number on the flow angles, dual vanes 
for both angle of attack and flank angle of attack should be 
used, making the flow field symmetrical vertically and lat- 
erally. Increasing the distance between the sets of angle- 
of-attack vanes and flank angle-of-attack vanes should re- 
duce the compression effect on the vanes, and eliminating 
the noseboom attachment collar will eliminate any compres- 
sion effect on the flank angle-of-attack vanes. 

On this program, problems were encountered with the 
flow angle potentiometers, which changed calibrations as 
they wore. The angle sensing elements of a syncro angle 
transducer do not rub against each other, and would not be 
susceptible to this problem. 

Since an accurate knowledge of the angular offset of the 
noseboom from the IRU is necessary for a successful cali- 
bration, these two systems should be rigidly attached to the 
same structure. The structure should have a location to mea- 
sure externally roll, pitch, and yaw angles on the ground. 
This would minimize the effects of maintenance and aircraft 
bending and torsion in flight on the misalignment angles. 

Concluding Remarks 

The research airdata system of an F-104 aircraft has been 
calibrated to measure winds aloft The wind measurement 
accuracy desired for this investigation required an airdata 
calibration more accurate than for most aircraft research 
programs. Mach number was calibrated using the tower 
fly-by and radar acceleration-deceleration techniques, with 
a flight-to-flight calibration repeatability of ±0.003 subsoni- 
cally and ±0.005 supersonically. Total temperature was cal- 
ibrated and found to have a recovery factor of 0.986 with a 
±0.009 scatter in the data. 

Flow angles were corrected for noseboom misalign- 
ment, time delays, aircraft angular rates, and noseboom 
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bending, then calibrated using trajectory reconstruction us- 
ing a multiple state linear Kalman filter (LKF). Nose- 
boom misalignment was corrected by an Euler rota- 
tion through its offset angles. The flow angles were 
corrected for aircraft angular rates without using small 
angle approximations. 

The trajectory reconstruction calibration method using 
the multiple-state LKF allows several different data sources 
to be blended, minimizing systematic errors in the data and 
yielding highly accurate trajectory information. The flight- 
to-flight calibration repeatability of angle of attack is ±0.2° 
subsonically and ±0.3° supersonically, and the flight-to- 
flight calibration repeatability of flank angle of attack is 
±0.15° subsonically and ±0.25° supersonically. These cal- 
ibration repeatabilities are considered to be near the attain- 
able accuracy limit with the configuration used. The sub- 
sonic flank angle-of-attack bias data had too much scatter 
to calibrate, probably because of the high sensitivity of the 
true flank angle-of-attack calculation to small errors in the 
meteorological wind estimate. 

Assuming negligible errors in the inertial velocities cal- 
culated by the inertial reference unit (IRU) or ground based 
radar, the airdata calibration repeatabilities attained should 
result in measured wind speed repeatabilities of 3 to 10 ft/sec 
(1 to 3 m/sec), approximately the same repeatability as the 
radar tracked Jimsphere balloon. 

Several recommendations from this investigation will be 
useful in the design and operation of future airdata systems 
for aircraft used to measure winds aloft. These include us- 
ing a noseboom with dual angle-of-attack and flank angle- 
of-attack vanes to reduce the sensitivity of upwash and side- 
wash on Mach number. The noseboom and IRU should be 
rigidly attached to the same structure to minimize geometric 
alignment variability. 
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Table 1. Noseboom misalignment. 


Flight Period 

<t>\ deg 

0\deg 

yi',deg 

Check 1 

-1.33 

-0.40 

0.53 

Check 2 

-1.25 

-0.23 

0.27 


Table 2. Sensor transformation distances. 

Distance from 
Parameter c.g., ft 

Distance from 
noseboom to 
accelerometer, ft 

x a 35.12 

y a -0.57 

xp 34.77 

zp 1.09 

m 

v_z 

-20.15 

0.2552 

-1.130 


Table 5. Flank angle-of-attack observation weights. 


Weighting 

Measurement 

North 

East 

Down 

Aircraft position 

10 

10 

10 

Inertial velocity 

1000 

1000 

1000 

Winds 

10 

10 

10 


Indicated 

Meteorological 

Mi 

Pressure 

1 

2 

2 


Tux 

Meteorological 


Temperature 

5 

1 



Table 6. Angle-of-attack dynamic equations weights. 




Weighting 


Measurement 

North 

East 

Down 

Aircraft position 

100 

100 

100 

Inertial velocity 

100 

100 

100 

Winds 

1 

10 

10,000 


Ambient 

Mach Position error 

Pressure 

1 

1 

1 

Temperature 

Ambient 

1 




Table 3. Time delays for IRU. 


Parameter 

Time delay, 
sec 

4> 

0.050 

9 

0.050 


0.110 

p 

0.050 

Q 

0.050 

r 

0.050 

True ground speed 

0.110 

True ground track angle 

0.110 

Flight path angle 

0.110 

Vb 

0.065 


Table 4. Angle-of-attack observation weights. 


Weighting 

Measurement 

Noth 

East 

Down 

Aircraft position 

10 

10 

10 

Inertial velocity 

100 

100 

100 

Winds 

10 

10 

0.001 


Indicated 

Meteorological 

Mi 

Pressure 

1 

1 

1 


Tux 

Meteorological 


Temperature 

2 

1 



Table 7. Flank angle-of-attack dynamic equations weights 




Weighting 


Measurement 

North 

East 

Down 

Aircraft position 

10 

10 

10 

Inertial velocity 

10 

10 

10 

Winds 

200 

200 

10,000 


Ambient 

Mach Position error 

Pressure 

1 

1 

1 

Temperature 

Ambient 

1 




Table 8. Total temperature recovery 
factor calibrations. 


Maneuver 

Recovery factor, k 

1 

0.9897 

2 

0.9953 

3 

0.9913 

4 

0.9766 

5 

0.9786 

Average k 

0.986 
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Fig. 1 The F- 104 aircraft. 



Fig. 2 Standard NACA airdata probe. 
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(b) Total and static port configurations. 
Fig. 2 Concluded. 



Fig. 3 Nondeiced, open element, total temperature sensor. 



Fig. 4 The IRU, AICS, PCM system, and signal condition 
ing installed in aircraft. 
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Fig. 5 Tower fly-by position error calibration method. 
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Fig. 6 Radar acceleration-deceleration position error calibration method. 
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Fig. 8 Flow angle correction steps. 
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Fig. 9 Flow angle definitions. 
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Fig. 10 Noseboom offset angles for check 1. 
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Fig. 19 Supersonic angle-of-sideslip calibration. 
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